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We present an updated and improved equation of state (which we call the NDL EoS) for use in 
neutron-star structure and supernova simulations. This EoS is based upon a framework originally 
developed by Bowers & Wilson, but there are numerous changes. Among them are: (1) a reformu- 
lation in the context of density functional theory; (2) the possibility of the formation of material 
with a net proton excess (Y e > 0.5); (3) an improved treatment of the nuclear statistical equilib- 
rium and the transition to heavy nuclei as the density approaches nuclear matter density; (4) an 
improved treatment of the effects of pions in the regime above nuclear matter density including the 
incorporation of all the known mesonic and baryonic states at high temperature; (5) the effects of 
3-body nuclear forces at high densities; and (6) the possibility of a first-order or crossover transition 
to a QCD chiral symmetry restoration and deconfinement phase at densities above nuclear matter 
density. This paper details the physics of, and constraints on, this new EoS and describes its im- 
plementation in numerical simulations. We show comparisons of this EoS with other equations of 
state commonly used in supernova collapse simulations. 



I. INTRODUCTION 

To describe the hydrodynamics of compact 
matter; be it in heavy-ion nuclear collisions, su- 
pernovae or neutron stars; an equation of state 
(EoS) is needed to relate the physics of the state 
variables. In supernovae the EoS determines 
the dynamics of the collapse and the outgo- 
ing shock, and determines whether the remnant 
ends up as a neutron star or a black hole. In a 
neutron star, it determines the maximum mass, 
mass-radius relationship, internal composition, 
cool-down time and dynamics of neutron star 
mergers. 

At present, only a few hadronic EoSs are com- 
monly employed that cover large enough ranges 
in density, temperature and electron fraction 
to be of use in core-collapse supernova simu- 
lations. These EoSs are usually in the form of 
a multi-dimensional table with three indepen- 
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dent variables (e.g. density, temperature and 
electron fraction). 

The two most commonly used equations of 
state in astrophysical simulations are the EoS 
of Lattimer & Swesty (LS91) pQ and that of 
Shen et. al. (Shen98) [3]. The former uti- 
lizes a non-relativistic parameterization of nu- 
clear interactions in which nuclei are treated 
as a compressible liquid drop including surface 
effects. The latter is based upon a Relativis- 
tic Mean Field (RMF) theory using the TM1 
parameter set in which nuclei are calculated 
in a Thomas-Fermi approximation. Recently, 
Shen et. al. [1] released a updates of the Shen98 
EoS table. The first update, EoS2, increased 
the number of temperature points as well as 
switching to a linear grid spacing in Y p . In the 
second update, EoS3, the effects of A hyper- 
ons were taken into account. It should be noted 
that several extensions to the Shen98 table have 
also been proposed, either by the implementa- 
tion of hyperons [5] or, of particular relevance 
to the present work, including a mixed phase 
transition to a quark gluon plasma [5J [7] . 

Over the last several years much progress has 
been made on other formulations of the super- 
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nova EoS. Several new EoS tables, each based 
upon RMF models, have been introduced [8j |9| . 
The new hadronic tables of Shen et. al. [5] are 
based upon a virial expansion and two differ- 
ent RMF interactions in the Hartree approxi- 
mation. The first table includes the NL3 [10] 
parameter set while the second table is param- 
eterized by the FSUgold parameters [TTj . The 
new EoS of Hempel et. al. [H] is described by an 
RMF in nuclear statistical equilibrium (NSE) 
for an ensemble of nuclei and interacting nucle- 
ons. Steiner et. al. [5] also constructed several 
new EoSs to match recent neutron star obser- 
vations. In these nucleonic matter was param- 
eterized with a new RMF model that treated 
nuclei and non-uniform matter with the statis- 
tical model of Hempel et. al. [T2"] . 

In this work we describe a new Notre Damc- 
Livermore (NDL) EoS that we make publicly 
available. This EoS evolves from the original 
Livermore formulation [HUH]. This NDL EoS, 
is consistent with known experimental nuclear 
matter constraints and recent mass and radius 
measurements of neutron stars. 

Below nuclear matter density, the conditions 
for NSE are imposed in the NDL EoS above a 
temperature of T « 0.5 MeV. Below this tem- 
perature the nuclear matter is a approximated 
by a nine element reaction network which must 
be evolved dynamically. Above this tempera- 
ture, the nuclear constituents are represented by 
free nucleons, alphas and a single "representa- 
tive" heavy nucleus. The high density phase of 
the EoS is treated with a parameterized Skyrme 
energy density functional that utilizes a modi- 
fied zero range 3-body interaction. The effects 
of pions and of the mesonic and baryonic reso- 
nances on the state variables at high densities 
are also included as well as the consequences 
of a phase transition to a QGP. For the NSE 
and high density regime, the EoS is provided in 
tabular form covering the necessary ranges in 
density, temperature and Y e needed for use in 
astrophysical simulations. Below nuclear satu- 
ration where NSE cannot be applied, numerical 
routines are available to the user [15] . 



II. THE NDL EQUATION OF STATE 

Depending upon the density and tempera- 
ture there are a variety of matter components 
that contribute significantly to the equation of 
state during various epochs of supernova col- 
lapse and the interiors of neutron stars. These 
include photons, electrons, positrons, neutrinos, 
mesons, all the 215 known mesonic and baryonic 
states [T5], free neutrons, protons, and atomic 
nuclei, and even the possibility of a crossover to 
a quark gluon plasma. At low density and high 
temperatures we assume a meson gas; consist- 
ing of thermally created, pair-produced mesons 
with zero chemical potential. In the high den- 
sity limit, but low temperatures, pions are con- 
strained by chemical equilibrium among the 
neutrons, protons and the other baryonic states. 
Baryons are assumed to have the same, non- 
zero chemical potential, such that baryon num- 
ber is conserved. The inclusion of the additional 
mesonic and baryonic states is yet another im- 
provement appearing in this updated EoS. 

Since the material is optically thick to pho- 
tons, one can include photons along with matter 
particles in the equation of state. The NDL EoS 
is divided into four regimes: 

1 . Baryons below nuclear matter density and 
not in NSE; 

2. Baryons below nuclear matter density and 
in NSE; 

3. Hadronic matter above saturation density 
including pions; and 

4. A first order or cross-over phase transition 
to quark gluon plasma. 

The electrons and positrons are approxi- 
mated as a uniform background and are treated 
as a non- interacting ideal Fcrmi-Dirac gas. Pho- 
tons are approximated as a black-body and thus 
are given by the usual Stefan-Boltzmann law. 
Neutrinos, however, are not necessarily confined 
and must be transported dynamically. In su- 
pernova simulations most matter except neutri- 
nos can be assumed to be in local thermody- 
namic equilibrium (one temperature in a zone) 
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but not necessarily in chemical equilibrium (i.e. 
the weak reactions have not necessarily equili- 
brated). The independent variables generally 
chosen for the equation of state are then the 
temperature T, the matter rest-mass density p, 
and the net charge per baryon Y e — n e /riB- The 
previous formulation required that Y e < 0.5, 
but we have removed that restriction in this new 
version. 



Baryons Below Saturation 
and not in NSE 



Below nuclear saturation density and above 
a temperature of T « 0.5 MeV we can as- 
sume that NSE is valid. Below this tempera- 
ture the isotopic abundances must be evolved 
dynamically. To achieve this, the nuclear con- 
stituents are approximated by a 9 element nu- 
clear burn network consisting of n, p, 4 He, 12 C, 
16 0, 20 Ne, 24 Mg, 28 Si, 56 Ni [17. The free en- 
ergy per baryon is taken to be the sum of con- 
tributions from an ideal gas Fg [Eq. ([!])] and a 
coulomb correction F c [Eq. JM]. The ideal gas 
contribution is simply, 
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The relevant variables here are the nuclear mass 
fraction X il the local baryon number density n, 
the temperature T, and the atomic mass num- 
ber Ai . The index i runs over the entire reaction 
network, and A is the thermal wavelength per 
baryon given by 



A = 



8tt 3 



e (27ttob) 
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[Note, that natural units (h = c = k = 1) 
have been adopted here and throughout this 
manuscript.] 

The Coulomb contribution is given by 



(3) 



From these relations the baryonic pressure and 
energy per unit mass can be calculated from the 
ideal gas thermodynamic relations. 
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B. Baryons Below Saturation and in NSE 

When NSE is valid, the baryonic nuclear ma- 
terial is approximated as consisting of a 4 com- 
ponent fluid of free protons, neutrons, alpha 
particles and an average heavy nucleus. This 
formulation is reasonably accurate and conve- 
nient in that it leads to fast analytic solutions 
for the NSE. One should exercise caution, how- 
ever, [T3] when considering detailed thermonu- 
clear burning or a precise value of Y e in NSE is 
desired. In such cases an extended NSE network 
should be employed. 

In the absence of weak interactions the neu- 
tron and proton mass fractions are constrained 
by charge conservation (i.e. constant electron 
fraction Y e ), 



YiXi — Y e 

i 

and baryon conservation, i.e. 

5> = 1 - 



(6) 



(7) 



The thermodynamic quantities are deter- 
mined from the Hclmholtz-frcc energy per 
baryon which is given as a sum of the various 
constituents, 



/ = /„ + fp + U + I A 



(8) 



The constituent free energies can be written an- 
alytically as 
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where the various terms in Eqs. ^ - (12 1 are 
defined as follows: 

(A) = 194.0(1 - r e ) 2 (l + X + 2X 2 + 3X 3 ), 

(13) 

is the density dependent mass of the average 
heavy nucleus, expanded in terms of the density 
parameter X, defined by 
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7.6 x 10 13 g cm 
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(14) 



X B is the free baryon mass fraction while X a 
and X A are the mass fractions of 4 He and the 
average heavy nucleus in obvious notation. The 
quantities Y p and Y n are the relative number 
fractions of free baryons in protons or neutrons, 
respectively. Thus, Y p + Y n = 1. The quan- 
tity Ya is the average ZIA for heavy nuclei. 
The quantity W in Eqs. |9| - (p~2]) is a weight- 



ing factor that interpolates between the low- 
density and high-density regimes. It is defined 
\yy W = (1 — p/pn) 2 . The transition from sub- 
nuclear to supra-nuclear density is expected to 
be continuous. The reason for this is that, as 
the density increases, the equilibrium continu- 
ously shifts to progressively heavier nuclei. 

When a relativistic Thomas-Fermi represen- 
tation of the electrons is evaluated at subnu- 
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clear density, the electron energy is lowered by 
more than ~1 MeV [13]. The electrostatic nu- 
clear energy also increases in magnitude. Simi- 
larly the transformation of nuclei from spheres 
to other more exotic shapes (e.g. pasta nu- 
clei, etc.) jT5J [TS] also lowers the energy of the 
medium by about 1 MeV. The net result is that 
the pressure and energy are smooth functions 
of density near the nuclear saturation density. 
Hence, the weighting factor is chosen to approx- 
imate this smooth transition. 

The normal 56 Fe ground state is taken as the 
zero of binding energy. This is unlike most other 
equations of state for which the zero point is 
chosen relative to dispersed free nucleons. The 
reason for the choice made here is that it avoids 
the numerical complication of negative internal 
energies in the hydrodynamic state variables at 
low temperature and density due to the bind- 
ing energy of nuclei. The energy per nucleon 
required to dissociate 56 Fe into free nucleons is 
e p0 = 8.37 MeV for protons, while for neutrons 
it is e„ = 9.15 MeV. 

The quantity p^ is the density at which nu- 
clear matter becomes a uniform sea of nucleons. 
This was found by fitting the saturation density 
of nuclear matter [i.e. Pm(p, T = 0, Y e ) = 0] as 
a function of p and Y e . The zero-temperature 
result was chosen to simplify the problem of 
making a smooth transition between the three 
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equation of state regimes. The result is 



taken to be 



p N = 2.66 x 10 
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(15) 



The quantities £ n and ( p in Eqs ^ and ( 10 1 
are a measure of the degeneracy of the free 
baryons. They are defined by 
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where the quantity B(pYiX B ) 2 ' 3 is the en- 
ergy per baryon of a zero-temperature, non- 
rclativistic ideal fermion gas and the constant 
B is 
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The dimensionless constant /3 appearing in 
Eqs. Q and (10) is determined such that the 
translational part of f p and /„ reduces to the 
correct non-degenerate limit (T — > oo, d —t 0). 
That is, 
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(19) 



where A is the thermal wavelength per baryon 
given m Eq. @. 

The function b(Y e ) in Eq. ( [12] ) is determined 
by the condition that the Coulomb contribution 
to the pressure at p = be canceled by the 
term proportional to b(Y e ). This requires, 
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The expression for the statistical weight of 
the heavy nucleus gA appearing in Eq. (12 1 is 




(21) 



where 



T s = (8 MeV) 1 
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(22) 



In Eq. (12) the constant S E = 120 MeV is de- 



rived for a symmetry energy of 30.4 MeV per 
nucleon (see below). The constant Yp e — 0.464 
is the fraction of protons in 56 Fe. 

The chemical potentials are found from the 
free energy as 
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(23) 
(24) 
(25) 
(26) 
(27) 



where p p , p n and p a are the chemical potentials 
of free protons, neutrons, and alpha particles. 
The quantities \i nA and p pA are the chemical 
potentials of neutrons and protons within heavy 
nuclei. These quantities are related by the Saha 
equation: 



2/^n ~\~ 2/ip Pa 
2p n A + 2p P A = Pa 
PnA ~ PpA = Pn ~ Pp = 



(28) 

(29) 
(30) 



In the original Livermore formulation |13l 114] , 
an analytical approximation was used to deter- 
mine the average heavy nucleus mass fraction, 
X A . In the current implementation, the three 
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chemical potential constraints combined with 
charge and baryon number conservation are 
solved self consistently to determine the matter 
composition. This leads to a 20% increase in the 
mass fraction of heavy nuclei when compared to 
the original approximation scheme |13[ 117) . 



C. Baryonic Matter Above Saturation 
Density 

Above nuclear matter density, the baryons 
are treated as a continuous fluid. In this regime, 
the free energy per nucleon is given in the form 

f = fx (n, Y p ) + f 2 (n, T) + 8.79 MeV , (31) 

where the addition of 8.79 MeV sets the zero for 
the free energy to be the ground state of 56 Fe. 
For an arbitrary proton fraction Y p and number 
density n the zero-temperature contribution to 
the free energy per nucleon is written as the sum 
of an isospin symmetric term and the symmetry 
energy: 

h («, Y p ) = - A (n, Y p = 0.5) + S (n, Y p ) . (32) 



Expanding S(n,Y p ) in terms of (1 — 2YJ,), and 
keeping only the leading contribution, the sym- 
metry energy can be written as 

S{n,Y p ) = {l-2Y p f So(n) , (33) 

where So can be identified as the symmetry en- 
ergy. 

Above saturation density we include both 2- 

(2*) (3) 

body (v\j ) and 3-body (v^jS) interactions in the 
many-nucleon system. The Hamiltonian of this 
system is thus given by 

£=£**+£4? + £ »S. (34) 

i i<j i<j<k 

where £$ is the one body contribution while Vij 
and Vijk are the 2 and 3-body interactions, re- 
spectively. In the density functional approach 
one can parameterize these interactions to de- 
scribe the ground-state properties of finite nu- 
clei and nuclear matter [2"0~H22"] . The micro- 
scopic interactions, such as meson exchange, are 
embedded in the parameters of the density de- 
pendent forces. 

Among the most widely used interactions are 
those of the Skyrme type forces. In this for- 
mulation the two-body potential is given in the 
form [23l : 
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to (l + xoPs) 6 (n - r 2 ) + iti (d (n - r 2 ) P + k 2 6 (r x - r 2 ) 
+ < 2 k 2 ■ 8 (tx - r 2 ) k + iW (a% + ct 2 ) • k' x <5 (n - r 2 ) k , 



(35) 



where P s is the spin exchange operator, r x and 
r 2 are the position vectors in the two-body po- 
tential, xo is the coefficient for the isospin ex- 
change operator, k and k' are the momentum 
and conjugate momentum operators, and Wq is 
the coefficient of the two-body spin orbit inter- 
action. 

We will discuss the Skyrme coefficients 
to> ti i ^3) and <J in the following sections. For 
this Skyrme potential the high density behavior 



can be dominated by a 3-body repulsive interac- 
tion. This term is taken to be a zero range force 
of the form v^z = fad (ri - r 2 ) 6 (r 2 — r 3 ). If 
the assumption is made that the medium is 
spin-saturated, which is valid for neutron star 
matter and nuclei |24j . the three-body term is 
equivalent to a density dependent two-body in- 
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teraction given by 



,(3) _ 



t 3 [l + P a )S (ri - r 2 ) n 



2 



In the present formulation we generalize this po- 
tential to a modified Skyrme interaction that 
replaces the linear dependence on the density 
by a power-law index <r. This modified Skyrme 
potential can then be written as |25j 



,0)' 



1 

= 6* 3 



l + P s U( ri -r 2 ) 



ri + r 2 



This modification has been introduced [2U] to 
increase the compressibility of nuclear matter 
at high densities. A value of a = 1/3 is a com- 
mon choice (26] [27]. However, in the present 
approach we choose to treat a as a free pa- 
rameter to be determined by constraining the 
third derivative of the energy per particle (e.g. 
skewness coefficient) from observed neutron- 
star properties [2"5] . 

The main advantage of the Skyrme density 
functional is that the variables that character- 
ize nuclear matter can be expressed as analytic 



functions. The isospin symmetric contribution 
is described by a Skyrme density functional with 
a modified three-body interaction term a. We 
use Tp to denote the kinetic energy of a particle 
at the Fermi surface 



T F = 
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(38) 



Then, calculating the expectation value of the 
Hamiltonian [Eq. (34)] in a Slater determinant 
and setting N=Z, the energy per nucleon for 
symmetric nuclear matter can be derived |23j . 



(37) as given in Eq. ( 39 1 



All quantities and coefficients for symmetric 
nuclear matter are obtained from Eq. ( |39| . The 

and 



pressure, is deduced from P 
is given in Eq. (40). Eq. (ill) gives the vol 




ume compressibility of symmetric nuclear mat- 
ter. This is calculated as the derivative of 
the pressure with respect to number density: 
# = 9(F) =18^+9n 2 

V on ) n 

skewness coefficient, Qo = 27n 



Finally, the 

3d 3 (E\ 
A) 



is de- 



duced from the third derivative of the free en 
ergy per nucleon jUj and is given in Eq. (142 
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(39) 
(40) 
(41) 
(42) 



These four equations completely describe the 
properties of symmetric nuclear matter. Values 
for the coefficients can be constrained by fix- 
ing the density of nuclear saturation, as well as 
imposing the observational constraint that the 



maximum mass of a neutron star must exceed 
1.97 ±0.04 M Q [28]. 

We employ a linearly increasing density de- 
pendent symmetry energy [301 131] . This be- 
havior is understood to arise from the use of a 
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vector coupled p-meson in relativistic field the- 
ory calculations. The form adopted here for the 



symmetry coefficient So is then 

72 



S (n) = ?y 



16 



1 + 4r/ 



(43) 



where rj = n/no is the saturation density pa- 
rameter with no the saturation number density 
(0.16 fm -3 ) [32]- Therefore, we find the zero- 
temperature contribution to the free energy to 
be 



3 3 1 

fi {n, Y p ) = -T F + -t Q n + y;t 3 n' 



<T + 1 
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D. Thermal Correction 

For the thermal contribution to the free en- 
ergy per particle we follow the approach de- 
scribed in Refs. [THE]. We assume a degener- 
ate gas of the baryonic states reviewed in [To] , 
as a function of temperature T and baryon den- 
sity n. Since the zero-temperature contribution 
to the free energy is already properly taken into 
account by the Skyrme and symmetry energy 
contributions, only the thermal portion needs 
to be added. We also assume that the baryonic 
states are in chemical equilibrium. The expres- 
sion for the thermal contribution is written: 

h{n, T) = a - a + - (w - w ) , (45) 
n 



where a and w are the finite temperature 
"chemical potential" and the grand potential 
density, respectively. The quantities, ao and ujq 
are the zero-temperature limits of the "chemical 
potential" and grand potential density and n is 
the local baryon number density. 



ao is constrained from the number density of 
baryons and is determined from the relation 



£^(«o-^) 3/2 , (46) 



while the zero-temperature limit of the grand 
potential density is 



w 



V — 

^ 2tt 2 



-rrii In 
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(47) 



In Eqs ( 46 ) and ( 47 ) m t is an effective particle 



mass deduced from fits to results from relativis- 



tic Bruckner Hartree-Fock theory 

mi 

1 + 0.27r? ' 



(48) 
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and the sum over i includes both nucleons and 
delta particles. 

The finite temperature a is then found from 
baryon number conservation 



E 



9i 
2tt 2 



(h(p, a) - h(p,-a))p 2 dp, 



(49) 

and the finite temperature grand potential den- 
sity is given as 

/*oo 4 

w = - ^2 ^2 / — ( h ( p > a ) + h & dp > 







where is given by 



e; = \ P 2 + mf- 



(50) 



(51) 



In Eqs. (49) and (50), h(p,a) is just the usual 



Fermi distribution function 

h{pi,a>i) 



exp[(e 2 - oci)/T] + 1 



(52) 



and (?i is the spin-isospin degeneracy factor. It 
should be noted that a and «o ar e only used 
to construct fa and do not correspond to a real 
chemical potential. The actual chemical poten- 
tials are found from derivatives of the total free 



energy [Eq. 31 with respect to density. They 
are given as 

-^ ! + "{l) T , Y ,- Y '(^) T „ <53) 



fl p = f + 71 

where Y n +Y p = 1 



E. Thermodynamic State Variables 

Once the thermal contribution to the free en- 
ergy is constructed the thermodynamic quan- 
tities can be calculated. Of particular interest 
are the total internal energy (E), the total pres- 
sure (P), the entropy per baryon (S) and the 



adiabatic index T. The total internal energy is 
calculated from the free energy: 



E = f-T 



dT 



E e + E. 
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(55) 



where E e is the electron energy determined by 
numerically integrating over Fermi-Dirac dis- 
tributions and E~ is the photon energy con- 
tribution determined from the usual Stefan- 
Boltzmann law. 

The pressure is calculated from 
P = n 2 (df /dn). It is important to note 
that the thermal contribution to the pressure 
is not the simple form of u>o — lu as one would 
expect from the usual application of the ther- 
modynamic potential, but is given by a slightly 
more complicated form 



d 

Ptherm = v - w + n— [uj - u> ] . (56) 



This is due to the fact that the effective mass 
[cf. Eq. (48)] is density dependent. If this de- 



pendence were removed we would recover the 
usual form of the pressure from the thermody- 
namic potential. 

We thus determine the total pressure from 
the free energy 



P 



dh 
dn 



+ P 



(lie 



(57) 



where again the electron and photon contribu- 
tions are determined from the pressures cal- 
culated previously. The entropy per baryon 
in units of Boltzmann's constant is given by 
a simple derivative S — — (df /dT), and the 
adiabatic index is given by the usual form 
r = [d InP/d ln(n)]. 

To investigate the physical quantities 
relevant for modeling nuclear matter in 
dense astrophysical environments we de- 
termine the coefficients in Eq. (39) by 



utilizing known constraints on nuclear matter 
saturation. Applying the saturation con- 
dition P(n = n ) = n 2 -^ (f ) | n=no = to 



Eqs. ( 39p2 ), one gets a system of four equa- 
tions in terms of the quantities t Q , (3ii + 5^), 
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*3, and a: 



En = 



-Tp n H — tnnn H t-inZ +1 + 

5 8 16 40 V 2 



3 /3^ 2 



2/3 



(3ti + 5*a) ^o /3 



2 3 1 , x „, , 1 f3n 2 

Po=0= -T Fo + -t n a + -ta (a + 1) n^ 1 + - ( — 



2\ 2/3 



(3ii + 5* 2 ) n 



9 3(¥) 



2/3 



5 Tfo + 16 i3Cr ^ + ^ n ° 



24 



27 



Qo = y r fb + Y6 <3Cr ( CT + x ) ( CT _ x ) n < 



<J+1 



(3ti + 5t 2 ) n 

2/3 



5/3 



5/3 



(3ti + 5< 2 ) n^ 3 



(58) 
(59) 

(60) 
(61) 



Solving Eqs. (58 1 - (61) for a then yields 
%T Fo - 2K + Qo - 45£o 



3if + 45£ 



27 rp 



(62) 



The usual approach is to choose a set of 
data, e.g. resonances, nuclear masses, charge 
radii, etc., to find a best set of parameters for 
a Skyrme model. From these Skyrme coeffi- 
cients, then, the quantities in Eqs. (58) - (61) 
are deduced, i.e. no, Tp , En, Kq, Qo- For our 
purposes, however, we choose a more empirical 
approach. That is, we adopt inferred values of 
no, Eg, Kq, Qo from the literature and use these 
to determine the Skyrme model parameters. We 
also demand that these parameters allow neu- 
tron star masses > 1.97 ± 0.04 M Q . In this way 
our EoS relates directly to inferred properties of 
nuclear matter in the literature and our EoS is 
easily adaptable to improved experimental and 
theoretical determinations. 

The saturation density no ~ 0.16 fm~ 3 and 
the binding energy per nucleon E = — 16 McV 
are reasonably well established [32] ■ The deter- 
mination of the compressibility parameter from 
experimental data on the giant monopole res- 
onance on finite nuclei, however, has been a 
long standing conundrum. On the one hand, 
the compressibility of nuclear matter can be 



Parameter 


value 


ref. 


n 


0.16 ± 0.01 fur 3 


m 


E 


-16 ± 1 MeV 


m 


K 


240 ± 10 McV 


m 


Qo 


-390 ± 90 MeV 


this work 



TABLE I. Adopted constraints on properties of nu- 
clear matter. 



determined |35j by fitting measured breathing- 
mode energies, using generalized Skyrme-type 
forces that include a density and momentum 
dependent term. Acceptable fits are in the 
range of K = 215 ± 15 MeV [35 . On the other 
hand a value of K = 231 ± 5 MeV, has been 
found 36 by using measured E0 distributions 
in 40 Ca, 90 Zr, 116 Sn, 144 Sm and 208 Pb based 
upon the calculations of [37]. Building a new 
class of Skyrme forces in Ref. [33] a value of 
Kq = 240 ± 10 MeV was found from these data. 

Given the unresolved discrepancy between 
GMR data and models [29], there is cur- 
rently a fairly large uncertainty in K . How- 
ever, for our purposes we adopt the me- 
dian value and uncertainty from Ref. [34) . i.e. 
A'o = 240 ± 10 McV as this is most appropri- 
ate for the Skyrme force approach employed 
here. Solving Eqs. (58) - (61) self consistently, 



we therefore determine the best range for the 
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FIG. 1. Mass vs. radius diagram constrain- 
ing the skewness coefficient Qo- The lower 
bound of Qo > -475 MeV is constrained by the 
M max > 1.97 ±0.04 M Q [28] observational limit. 



nuclear compressibility consistent with the re- 
sults of [33]. 

There is even more uncertainty in the skew- 
ness parameter Qq. In Ref. [35] . breathing 
mode data were used to find a weak inverse 
correlation between Qo and K . In that work 
however, they could only deduce a very broad 
range for Q = -700 ± 500 MeV. We find an 
upper bound on the skewness coefficient by 
using the stiffest compressibility in our range 
(Kq = 250 MeV) in Eqs. d58b - (|6lJ. This gives 



us an upper bound on the skewness coefficient 
of Qq = -300 MeV. The lower bound on Q 
is found similarly, using the smallest compress- 
ibility in our range (Kq — 230 MeV). It is fur- 
ther constrained that the maximum mass of 
the neutron star is above the 1.97 ± 0.04 M 
observation. From Fig. [TJ the lowest skew- 
ness coefficient that meets this constraint is 
Qo = -387 ± 87 MeV. This is also consistent 
within the range given in Ref. [31]. The fidu- 
cial NDL EoS is constructed using the median 
values of the constraints given from Table [I] 



NDL EoS are then found to be 

t = -1718 MeV fm 3 
(3ti + 5*2 ) = -102 MeV fm 5 
t 3 = 13226 MeV fm 3<T+3 
cr = 0.369 . 

Note that our deduced value of the 3-body in- 
dex, a, is close to the commonly employed value 
of 1/3 [261 [27]. 

The density dependence of the symmetry en- 
ergy beyond saturation is highly uncertain. For 
many Skyrme models the symmetry energy ei- 
ther saturates at high densities, or in the worst 
case becomes negative. This results in a nega- 
tive pressure deep inside the neutron star core. 
For this work, we choose a fairly stiff symme- 
try energy, Eq. (43 1. That is, we implemented 



a linearly increasing function of density to re- 
move the issues inherent to many Skyrme pa- 
rameter sets. The symmetry energy at sat- 
uration is known [38] to lie within the range 
So = 26 — 34 MeV and is determined by the dif- 
ference between the energy per particle for pure 
neutron matter and that of symmetric mat- 
ter at T = MeV. For all relevant parameter 
sets chosen the NDL EoS symmetry energy at 
saturation is found to be So = 30.4 MeV (see 
Fig.g). 



III. PIONS IN THE NUCLEAR 
ENVIRONMENT 



The impact on the hadronic EoS from the 
lightest mesons (i.e. the pions) has been con- 
strained |39j from a comparison between rel- 
ativistic heavy-ion collisions and one-fluid nu- 
clear collisions. The formation and evolution 
of the pions was computed in the context of 
Landau-Migdal theory [33 to determine the 
pion effective energy and momentum. In this 
approach the pion energy is given by a disper- 
sion relation [33] 



The Skyrme coefficients for the fiducial 
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FIG. 2. The energy per baryon for pure neutron 
matter and symmetric nuclear matter. The sym- 
metry energy is calculated as the difference be- 
tween these functions at saturation density. For the 
NDL EoS a symmetry energy of So = 30.4 MeV was 
found. 



where fa„ is the pion "effective mass" defined 
to be 



i v = m„ \J\ + n(e 7r ,p Tr ,n) 



(64) 



Following [ST] and [ID] the polarization pa- 
rameter II can be written, 

where the denominator is the Ericson-Ericson- 
Lorentz-Lorenz correction |41j . The quantity 
A = exp(— p 2 /b 2 ) with b — Im^, is a cutoff that 
ensures that the dispersion relation [Eq. (63)] 



asymptotically approaches the high momentum 
limit, 



e (p w >■ co, n) = \Jm\ +pl~m N . (66) 



Following [UJ we take the polarizability to be 



(67) 



where a — 1.13 /m 2 . This form for the polar- 
izability ensures that the effective pion mass is 
always less than or equal to the vacuum rest 
mass m^. 

A key quantity in the above expressions is 
the Landau parameter g' . This is an effective 
nucleon-nucleon coupling strength. To ensure 
consistency with observed Gamow- Teller tran- 
sition energies a constant value of g' = 0.6 was 
used in [3D]. However, in [3DJ, Monte-Carlo 
techniques were used to statistically average a 
momentum dependent g' with particle distribu- 
tion functions. It was found that g' varies lin- 
early with density and is approximately given 
by 



9 = 9i + 92V 



(68) 



A value of gi = 0.5 was chosen to be consistent 
with known Gamow- Teller transitions. A value 
for <?2 was then obtained [39] by optimizing fits 
to a range of pion multiplicity measurements 
obtained at the Bevlac [32]. These data were 
best fit for a value of gi =0.06. 

The pions are assumed to be in chemical equi- 
librium with the surrounding nuclear matter. 
We consider the pion-nucleon reactions: 



p -H- n + Tt~ 



n o p + 7r 



(69) 



This leads to the following relations among the 
chemical potentials for neutrons, protons, and 
pions 



fj, p — [i n + fj,^+ , [i n — fi p + fi^- 



(70) 



These equilibrium conditions let us express the 
pion chemical potentials in terms of the neutron 
and proton chemical potentials: jl = \x n — [i v — 
/i^- = . Using the definitions of /i„ and 

Hp from Eqs. (53 1 - ( 54 ), the expressions for the 
pion chemical potentials are found to be 



4ff (n,y p ) 
(1 - 2Y p ) 



(71) 



Where S{n, Y p ) is the nuclear symmetry energy 
from Eq. ([43| . 



For a given temperature (T) and number den- 
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sity (n) the pion number densities are given by 
the standard Bose-Einstein integrals 



dp 



2ir 2 e^-^)/ T - 1 



(72) 



where i sums over {ir + , ir , 7r }, and is given 
by Eq. ( 63 ). Note that the ir° chemical potential 



is taken to be zero, since these particles can be 
created or destroyed without charge constraint. 

The charge fraction per baryon for the 
charged pions is defined as Y n - = n^- jn. From 
Eq. ( 71 ) we can calculate the pion number den- 



sities from the pion chemical potentials. Then, 
electric charge conservation gives, 



Y — Y — Y 



(73) 



Thus, we can solve Eq. ( 73 ) for the unknown 
quantity Y p . 

Once Y p is determined, the pionic energy den- 
sities and partial pressures can be calculated 
from 



E, = 
and 
Pi = 



V 



2ir 2 exp [(e„ -m)/T\-l ' 



p 2 (l/3)p(dejdp) 
2tt 2 exp [( Cff - Mi) /T] - 1 



(74) 



(75) 



Note that in the high temperature, low-density 
regime we add all baryonic and mesonic reso- 
nances. In this limit the pionic mass approaches 
the bare pion mass. Hence, we also trial all 
mesonic and baryonic states using bar masses. 



IV. QCD PHASE TRANSITION 

It is generally expected [13] that for suffi- 
ciently high densities and/or temperature, a 
transition from hadronic matter to quark-gluon 
plasma (QGP) can occur. Recent progress [44] 
in lattice gauge theory (LGT) has shed new 
light on the transition to a QGP in the 
low baryochemical potential, high-temperature 
limit. It is now believed that at high tem- 



perature and low density a deconfinement and 
chiral symmetry restoration occur simultane- 
ously at the crossover boundary. In particu- 
lar, at low density and high temperature, it 
has been found [H] that the order param- 
eters for deconfinement and chiral symmetry 
restoration changes abruptly for temperatures 
of T = 145 - 170 MeV g5] HB]. However, nei- 
ther order parameter exhibits the characteristic 
change expected from a 1st order phase tran- 
sition. An analysis of many [57J H5J thermo- 
dynamic observables confirms that the transi- 
tion from a hadron phase to a high temperature 
QGP is a smooth crossover. 

At low density the hadron phase can be ap- 
proximated as a pion-nucleon gas, while the 
QGP phase can be approximated as a non- 
interacting relativistic gas of quarks and glu- 
ons |49j . Equating the pressures in the hadronic 
and QGP phases, the critical temperature T c 
for the low density transition can be approxi- 
mated [Ml as: 



T c ~ {g q 



91, 



,-1/4 



1/4 



(76) 



Where the statistical weight g q for a low-density 
high-temperature QGP gas with three relativis- 
tic quarks is g q « 51.25, while gh ~ 17.25 was 
found for the hadronic phase by summing over 
all known meson data. 

Adopting the lattice gauge theory results [44] 
that 145 < T c < 170 MeV, then implies [49] 
that a reasonable range for the QCD vacuum 
energy is 165 < B 1 / 4 < 240 MeV. This pro- 
vides an initial range for the QCD vacuum en- 
ergy. We will further constrain this parameter 
by requiring that the maximum mass of a neu- 
tron star exceed 1.97 ± 0.04 M [25] . 

Another parameter that impacts the thermo- 
dynamic properties of the system is the strong 
coupling constant a s . For this manuscript we 
adopt a value of a s = 0.33 as this is a rep- 
resentative value for the energy regime under 
consideration [T5] . 

A transition to a QGP phase during the 
collapse can have a significant impact on the 
dynamics and evolution of the nascent proto- 
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neutron star. In [5U] it was first shown that a 
first order phase transition to a deconfincd QGP 
phase resulted in the formation of two distinct 
but quickly coalescing shock waves. More re- 
cently, it has been shown [5T] that if the transi- 
tion is first order, but global conservation laws 
are invoked, then the two shock waves can be 
time separated by as much as ~ 150 ms. Neu- 
trino light curves showing such temporally sepa- 
rated spikes might even be resolvable in modern 
terrestrial neutrino detectors [? ]. 

The observation of a 1.97 ± 0.04 M neutron 
star, however, constrains the possibility of a 
first order phase transition to a quark gluon 
plasma taking place inside the interiors of stable 
cold neutron stars [25]. Nevertheless, for initial 
stellar masses beyond > 20 Mq every phase of 
matter must be traversed during the formation 
of stellar mass black holes. Hence at the very 
least, this transition to QGP may have an im- 
pact |52j on the neutrino signals during black 
hole formation as well as its possible impact on 
core-collapse supernovae. 

A. The Quark Model 

For the description of quark matter we uti- 
lize a bag model with 2-loop corrections, and 
construct the EoS from a phase-space integral 
representation over scattering amplitudes. We 



allow for the possibility of a coexistence mixed 
phase in a first order transition, or a simple di- 
rect cross over transition. In the hadronic phase 
the thermodynamic state variables, are calcu- 
lated from the Helmholtz free energy F(T, V, N) 
as described in the previous sections. However, 
it is convenient to compute the QGP in terms 
of the grand potential, S1(T, V, ju). Both descrip- 
tions are equivalent and are related by a Legen- 
dre transform: SI = F — J2i Mi-^i- 

The grand potential for the quark-gluon 
plasma takes the form: 

n = J2(n q0 + n q2 ) + n g o + n g2 + BV. (77) 

i 

Where go and go denote the th -order bag model 
thermodynamic potentials for quarks and glu- 
ons, respectively, while q 2 and g 2 denote the 
2-loop corrections. In most calculations suffi- 
cient accuracy is obtained by using fixed cur- 
rent algebra masses (e.g. m u ~ m d ~ GeV, 
m s ~ 0.1 — 0.3 GeV). For this work we chose 
the strange quark mass to be m s = 150 MeV 
and a bag constant B 1 / 4 = 165 - 240 MeV. The 
quark contribution to the thermodynamic po- 
tential is given [33] in terms of a sum of the 
ideal gas contribution plus a two loop correc- 
tion from phase-space integrals over Feynman 
amplitudes [53]: 



2N C T 



d 3 p 



Sl q2 =« S 7T {N 2 C - 1 



„ (2^ 2 ) 

i r°° d 3 P Ni( P ) 



In ( 1 + e-W'-rt) + In ( 1 + e ~^ + ^ 

1 



3j a {2^)E i {p)'} (27T*) {2^)E l {p)E l {p') W(pW) + 2] 
N+ (p)N+ {p') + Nr ( p ) N r (p>) N+ (p)N+ (j,>) + Nr ( p ) N r (p') 



{E^-E^p')) 2 -{p-p'f 



{E^p) - E^p')) 2 ~{p-p>) 2 



(78) 



(79) 



where the denote Fermi-Dirac distributions: 

1 



eiS^WwO + l 



(80) 



r 



The one- and two-loop gluon and ghost con- 
tributions to the thermodynamic potentials can 
be evaluated in a similar fashion to that of the 
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quarks. 



/>oo t3 



45 



iV 2 



IT 4 







92 



3G 



a s 7V c (iV c 2 - 1) T 4 



(81) 



(82) 



For the massless quarks, Eqs. ( 78||79 ) are eas- 
ily evaluated to give 



2tt ! 



z _ (N?-l)a s /5tt 2 4 , 2 M * 
S ^ 2 ~ 8^ I IT* + 



(84) 



For the massive strange quark Eq. (78) can be 



easily integrated. Eq. (79 1, however, cannot be 



integrated numerically, due to the divergences 
inherent in it. We therefore, approximate |43) 
the two loop strange quark contribution with 
the zero mass limit. This may over estimate 
the contribution due to a finite strong coupling 
constant, but given that the quark mass is rela- 
tively small compared to its chemical potential, 
this is a reasonable approximation. 



B. Conservation Constraints 

The neutronized matter deep inside the core 
of a collapsing star consists of a multicompo- 
nent system constrained by the conditions of 
both charge and baryon number conservation. 
The pressure varies as a function of density 
for a first order transition producing a mixed 
phase of material. In fact, all thermodynamic 
quantities vary in proportion to the volume 
fraction [\ = V Q /(V Q + V H )] throughout the 
mixed phase regime. 

For the description of a first order phase tran- 
sition we utilize a Gibbs construction. In this 
case the two phases are in equilibrium when the 
chemical potentials, temperatures and the pres- 



sures are equal. For the description of the phase 
transition from hadrons to quarks this construc- 
tion can be written 



Ms 
Tq 



(85) 
(86) 
(87) 
(88) 



P H (T, Y e , {m? }) = PO (T, Y e , { Ml 9}) ,(89) 

where {pif} = {n n ,fJ. p , Me, AV} and {fif} = 
{/i tt ,W,/j s ,^,/i„}. 

The Gibbs construction ensures that a uni- 
form background of photons and leptons exists 
within the differing phases. Therefore, the con- 
tribution from the photon, neutrino, electron, 
and other lepton pressures cancel out in phase 
equilibrium. Also from this, we find that the 
two conserved quantities vary linearly in pro- 
portion to the degree of completion of the phase 
transition i.e. 

n B Y e =(l~ X )n%Y c H + X n%Y? (90) 
n B ={l- X )n% + X n% (91) 

where we have defined Y C H = Y p + Y„+ — Y n - 
and n®Y® = 1/3 (2n u — n^ — n s ). The inter- 
nal energy and entropy densities likewise vary 
in proportion to the degree of phase transition 
completion 



= (l-x)e H + X e Q (92) 
= (l-x)s H + xs Q - (93) 



V. RESULTS AND COMPARISONS 

In this section we compare properties of the 
new NDL EoS with the two most commonly 
employed equations of state used in astrophys- 
ical collapse simulations as well as the original 
EoS of Bowers & Wilson. Fig. [3] shows the to- 
tal internal energy as a function of local proper 
baryon density. We compare the Lattimer & 
Swesty EoS pQ at two differing compressibilities 
(K = 180 MeV and K = 220 MeV) with the 
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FIG. 3. The energy per particle as a function of 
density comparing the Lattimer & Swesty EoS with 
compressibilities K = 180 MeV and K = 220 MeV 
with the fiducial NDL EoS. 



FIG. 4. The pressure as a function of density com- 
paring two EoSs from Lattimer & Swesty pQ, the 
Shen EoS [J and the fiducial NDL EoS of the 
present work. 



fiducial NDL EoS at a fixed electron fraction 
and temperature of Y e — 0.3 and T = 10 MeV. 
A steep rise in the energy per baryon at high 
densities occurs for larger values of the com- 
pressibility as expected. 

Similarly, Fig. [4] depicts the pressure vs. den- 
sity for the Shen EoS [I] , the Lattimer & Swesty 
EoS [I] and the NDL EoS. The Shen EoS con- 
sistently leads to higher pressure. This results 
from the use of the TM1 parameter set that con- 
tains relatively high values for both the symme- 
try energy at saturation and the nuclear com- 
pressibility typical of RMF approaches. 



Pion Effects on the EoS 



The solution to the pion dispersion relation, 



Eq. (63), does not produce conditions within 



the supernova core to generate a pion conden- 
sate. The pions considered here are thermal pi- 
onic excitations calculated from the pion prop- 
agator in Eq. (65). In very hot and dense nu- 



clear matter the number density of pionic ex- 
citations is greatly enhanced by the irN A cou- 
pling |40j . Hence, it becomes energetically fa- 
vorable to form pions in the nuclear fluid when 



the chemical balance shifts from electrons to 
negative pions. This allows the charge states 
to equilibrate with these newly formed bosons. 
Since the pions are assumed to be in chemical 
equilibrium with the surrounding nuclear fluid, 
this has a profound effect on the proton fraction 
within the medium, particularly for low electron 
fractions, (see Fig. [5]). 

The pion charge fraction as a function of 
baryon number density is shown in Fig. [6j For 
a low fixed Y e the charge fraction of negative pi- 
ons can actually become greater than the elec- 
tron fraction, and the negative pions essentially 
replace the electrons in equilibrating the charge. 

From the solution to the pion chemical po- 
tentials [Eq. ( 71 1] , one finds that as the den- 
sity increases, negative pions are created due 
to the chemical potential constraints and the 
dispersion relation [Eq. (63)]. At the same 



time the number density of positively charged 
pions remains negligible due to the fact that 
it has a negative chemical potential. Due to 
its dependence on both the symmetry energy 
[Eq. (43)] and the isospin asymmetry parame- 
ter I = (1 — 2Y p ), the pion chemical potential 



[Eq. ( [71] )] increases linearly with respect to den- 
sity but decreases linearly with respect to Y p . 
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here, pions would not exist in the ground state 
configuration of a cold neutron star. As the 
temperature approaches zero the pionic effects 
diminish, until only the nucleon EoS contributes 
to the neutron star structure. 

In the hot dense medium of supernovae, how- 
ever, these pions tend to soften the hadronic 
EoS since they relieve some of the degeneracy 
pressure due to the electrons. We have found 
that the reduction in pressure is relatively in- 
sensitive to the temperature of the medium and 
is lowered by approximately 10% for all repre- 
sentative temperatures found in the supernova 
environment, as shown in Fig. [7] 



FIG. 5. The proton fraction above nuclear satura- 
tion showing the effects of pions in the hot dense su- 
pernova environment. For small electron fractions 
more pions are created due to the dependence of 
the chemical potential on the isospin asymmetry 
parameter /. 
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FIG. 7. Pressure versus baryon number density 
showing a 10% reduction in the pressure at high 
densities due to the presence of pions. 
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FIG. 6. Pion charge fraction versus density at a 
temperature T = 10 MeV. At a density of about 
0.65 fm -3 , the pion charge fraction exceeds 
the electron fraction of the medium. 



This will affect SN core collapse models since 
it allows collapse to higher densities and tem- 
peratures without violating the requirement 
that the maximum neutron star mass exceed 
1.97 ± 0.04 Mq for cold neutron stars. 



B. Hadron QGP Mixed Phase 



Therefore, for high electron fractions the pion 
chemical potential will remain small and charge 
equilibrium can be maintained solely among the 
electrons and protons. 

It should be noted, however, that as treated 



The constraint of global charge neutrality ex- 
ploits the isospin restoring force experienced by 
the confined hadronic matter phase. This por- 
tion of the mixed phase becomes more isospin 
symmetric than the pure phase because charge 
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is transferred from the quark phase in equilib- 
rium with it. 

Fig.[8]shows the charge fractions of the mixed 
phase and hadronic phases. From this we 
see that the internal mixed phase region of a 
hot, proto- neutron star contains a positively 
charged region of nuclear matter and negatively 
charged regions of quark matter until a density 
of no ~ 1.0 fm -3 . The presence of the isospin 
restoring force causes the thermal pionic con- 
tribution to the state variables to be negligible. 
This is due to the dependence of the pion chem- 
ical potential on the isospin asymmetry param- 
eter (1 — 2Y p ). As the hadronic phase becomes 
more isospin symmetric, the pion chemical po- 
tential remains small compared to its effective 
mass. 



tonically increasing function of density. 

Fig. [9] shows pressure versus density for vari- 
ous values of Y e , through the mixed phase region 
into a phase of pure QGP. One of the features 
shown is that as the density increases through 
the mixed region the slope of the pressure de- 
creases slightly. This becomes more evident 
when the adiabatic index, T, is analyzed as a 
function of density as shown in Fig. |10| Here, we 
find that the EoS softens abruptly upon enter- 
ing the mixed phase due to the fact that increas- 
ing density leads to more QGP rather than an 
increase in pressure. If this occurs while form- 
ing a proto-neutron star, its evolution will be 
affected as T falls below the stability point of 
T < 4/3 for n ~ 1.1 far 3 . 
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FIG. 8. Charge fractions of the mixed quark and 
hadronic phase. Due to the redistribution of charge, 
the hadronic phase becomes isospin symmetric even 
exceeding Y p > 0.5. This has the effect of lowering 
the symmetry energy and thus reducing the pres- 
sure in the hadronic phase. 

Since stars contain two conserved quantities, 
electric charge and baryon number, the coexis- 
tence region cannot be treated as a single sub- 
stance, but must be evolved as a complex mul- 
ticomponent fluid. It is common in Nature to 
have global conservation laws and not necessar- 
ily locally conserved quantities. Hence, within 
the Gibbs construction the pressure is a mono- 



FIG. 9. Pressure as a function of baryon number 
density through the mixed phase transition. The 
EoS softens significantly upon entering the mixed 
phase due to the larger number of available degrees 
of freedom. 

The collapse simulations of [50? ] show that 
as T falls below 4/3, a secondary core collapse 
ensues. The matter sharply stiffens upon enter- 
ing the pure quark phase at n ~ 1.4 fm -3 , and 
a secondary shock wave is generated. As this 
shock catches up to the initially stalled accre- 
tion shock, a more robust explosion ensues. 

Another feature seen in Fig. [9] is the Y e de- 
pendence of the onset density of the mixed 
phase. Fig. [TT] shows a phase diagram indi- 
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FIG. 10. F as a function of baryon number density 
showing the softening of the EoS as it enters the 
mixed phase regime. The EoS promptly stiffens as 
it exits the mixed phase into the pure quark matter 
phase due to losing the extra degrees of freedom 
supplied by the nucleons. 
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FIG. 11. Density-temperature phase diagram show- 
ing the density range of the mixed phase coexistence 
region for two values of Y e (dashed line for Y e =0.1 
and dash-dotted line for Y e = 0.3). The onset of the 
mixed phase is indicated by the set of curves on the 
left, while the curves on the right show the comple- 
tion of the mixed phase. For low temperatures it is 
seen that the onset density is highly Y e dependent. 
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FIG. 12. Neutron star mass-radius relation for var- 
ious values of the bag constant B 1 / 4 . We find that 
a first order phase transition is consistent with the 
maximum mass neutron star measurement for our 
adopted value of B 1/4 = 180 MeV. 



as a function of density for two values of Y e . 
For higher temperatures the onset happens at 
lower densities as would be expected. However, 
for high electron fractions (Y e ~ 0.3) such as 
those that can be found deep inside the cores 
of a proto-neutron star, the transition density 
remains quite high n c ~ 0.6 fm~ 3 . It is also 
of note that the coexistence region slightly de- 
creases as the electron fraction is increased. 



C. Neutron Stars with QGP Interiors 

As stated previously the observation [5S] 
of a 1.97 ± 0.04 M® neutron star has ruled 
out many exotic EoSs including many Hyper- 
onic models [3S]. However, using the range 
of bag constants determined by Eq. ( 76 1 we 
find that a first order phase transition to a 
QGP is consistent with the high maximum neu- 
tron star mass constraint |28| for our fiducial 



NDL EoS. In Fig. 12 we show that a bag con- 
stant B 1 / 4 > 190 MeV is required to satisfy the 
maximum neutron star mass constraint. This 
imposes a low baryon density transition tem- 
perature of T c > 150 MeV which is slightly 
below the current range of crossover temper- 
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atures determined from LGT jSJ. Hence all 
allowed values of the Bag constant inferred 
from LGT are consistent with the neutron 
star mass constraint. For our purpose we 
will adopt B 1 / 4 = 190 MeV (corresponding to 
T c ~ 150 MeV). We note, however, that the 
maximum mass is not above the observational 
limit if you simply ignore the 2-loop corrections 
(while keeping the same choice for B 1 / 4 ). 

Fig. [13] compares the neutron star mass ra- 



dius relation for the NDL EoS for: 1) a hadronic 
EoS (solid line); 2) a first order QCD transition 
with B 1 / 4 = 190 MeV (dot-dot dashed line); 
and 3) a simple QCD cross over transition (dot- 
ted line). Also, shown for comparison are re- 
sults from the LS220 (dashed line), LS180 (dot- 
dash line) Shen EoS (long-dashed line) and the 
original Bowers & Wilson EoS (dash-dash dot- 
ted line). Note, that all three versions of the 
NDL EoS easily accommodate a maximum neu- 
tron star mass > 1.97 ± 0.04 M w . 
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■ — ■ — LSI8OE0S 
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FIG. 13. Mass-radius relation for the Shen (dash-dash-dashed line), Lattimer & Swesty 180 & 220 (dashed 
line and dot-dash-dot-dash), Bowers & Wilson EoS (dash-dash dotted line), the NDL EoS with (solid line) 
and without (solid line) a mixed phase transition to quark gluon plasma as well as a simple crossover 
transition (dotted line) to a QGP. Note that all of these curves satisfy the 1.97 ±0.04 Mq astrophysical 
constraint except for the Bowers & Wilson and LS180. 
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VI. CONCLUSION 



We have discussed a much updated and im- 
proved equation of state based upon the orig- 
inal Livermore framework [13j [14]. We have 
shown that it is complementary to the most 
frequently employed equations of state for core- 
collapse supernovae due to Shen et. al. [2H3, 
and Lattimer & Swesty pQ. This NDL EoS is 
consistent with the known constraints of sym- 
metric nuclear matter and observed proper- 
ties of neutron stars and pulsars, whereas the 
previous version |T3l [14] was not. We found 
that consistently applying the constraints on 
symmetric nuclear matter, combined with the 
observation of a 1.97 ±0.04 Mq neutron-star, 
places a stronger limitation on the Skewness 
coefficient Qo than is available in the litera- 
ture. A first order phase transition to a QGP 
phase was also discussed in the context of a 
Gibbs construction. Applying the constraints 
from LGT for the range of low-baryon-density 
crossover temperatures, we were able to match 
the known constraints of the current maximum 



neutron star mass measurement for a bag con- 
stant B 1 / 4 > 140 MeV. On the other hand, if 
there is a cross over QCD transition the neu- 
tron star mass constraint can be easily accom- 
modated for any value of the bag constant. 

This confirms that a core collapse explosion 
paradigm, including a transition to quark gluon 
plasma may impact the neutrino light curve, 
shock dynamics, and heavy element nucleosyn- 
thesis via the ^-process and vp process both in 
supernovae and/or black hole formation. The 
consequences of this new NDL EoS for the dy- 
namics of core collapse supernovae, along with 
its impact on nucleosynthesis, will be explored 
in forthcoming manuscripts. 
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